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Abstract 

Thermodynamic properties of a fermionic system close to the unitarity Hmit, where the 2-body scattering 

\Q '• length a approaches ±00 is studied at high temperature Boltzmann regime. For dilute systems the virial 

O ! expansion coefficients in the Boltzmann regime are expected, from general arguments, to be universal. A 

pg I model independent finite temperature T calculation of the third virial coefficient b^{T) is presented. At the 

^ ' unitarity limit, b^ « —0.48 is an universal number. The energy density up to the third virial expansion is de- 

^' rived. These results are relevant for dilute neutron matter and could be tested in current atomic experiments 

on dilute fermionic gases near the Feshbach resonance. 
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Recent experiments in atomic traps near the Feshbach resonance 111 S |3l |4l |5|] have opened 
opportunities to study properties of finite density systems that have not been explored before. 
Using an external magnetic field, it is possible to tune the 2-body 5-wave scattering length a 
essentially at will: it can be made arbitrarily large (a = ±0°) compared to the range of the 2-body 
interaction R. Several groups have taken measurements near the Feshbach resonance in ^Li and 
^^Ka fermionic atomic gases llll|2l|3|,|4|,|5|]. These experiments are important checks of theoretical 
tools, both new and old, in previously untested territories and hold the prospect of revealing many 
new phenomenon. Equally important, what we learn from these atomic studies have universal 
application in many other subfields of physics. 

At densities n that are dilute compared to the range of the interactions {nF? <^ 1), the physics 
should be insensitive to the details of the interaction at the short distance scale R. Therefore, near 
the Feshbach resonance, as \a\ ^00, where there aren't any relevant scale left in the interaction, 
its properties are expected to be universal and not just applicable to only atomic systems. Many 
physical systems in nature are close to this universal limit. In "^He atomic gases, the scattering 
length a ~ lOOA is much larger than the range of the interaction 7? ~ 5A. To take an example from 
nuclear physics, the neutron-neutron 5- wave scattering length a„„ ~ — 19 fm is much larger than 
the range of the interaction set by the pion mass R ~ h/mTi ~ 1.4 fm. The universal properties 
learned in atomic experiments are applicable to other problems in the same universality class that 
are otherwise not directly accessible such as neutrino physics in supernova which depends on 
properties of dilute nuclear matter with a large resonating scattering length |6]. At the typical 
temperatures T ~ 4 MeV ILZllsl] and densities n ~ 10^^ fm^ ||6|] of the neutrinosphere, it is possible 
to calculate the equation of state to a given order in nh? expansion, where X = -sJl'KJiMT) is the 
non-relativistic thermal wavelength of a fermion with mass M. We calculate the equation of state 
up to the third order in nK^ and it model independently describes dilute fermionic systems at finite 
temperature with a large scattering length even if the microscopic physics is difficult to calculate 
or poorly understood. 

In this paper, we study the properties of a spin- 5 non-relativistic fermion at finite temperature 
T such that the thermal wavelength A, = ^/Iti/JmT) ^ R and we are not sensitive to the short 
distance scale R. However, the temperature is large enough such that the thermal wavelength 
is small compared to the inter-atomic distance nX^ <^ 1 and scattering length A, < a. Near the 
Feshbach resonance with \a\ — ;► 00, the hierarchy of momentum scales would be then \a\ ^ n~^'^ ^ 
X^ R. The calculation is more general and applicable for nX^ <^l,X^ R. 

The calculation is organized as an expansion in nX^ <^ 1, the Boltzmann regime, with factors 
of \a\y^MT /{In) summed to all orders in perturbation. In the Boltzmann regime the pressure P 
for a 2-component fermi system can be written in terms of the so called virial expansion ^: 

^ = ^[blZ + b2Z^ + hz' + ...]. (1) 

where z = exp{^/T) is the fugacity for the system with chemical potential /j. In this expansion, all 
the density dependence is in the fugacity z- This is a valid expansion for z ^ 1 or equivalently for 
nh? <^ 1 as shown later in Eq.fTSl 

The dimensionless virial coefficients bn depend on the vacuum interaction and the thermal 
momentum \JMT . For a dilute system with thermal wavelength X much larger than the range R, 
the interaction depends only on the scattering length a. Thus, as a ^ ±00 at the unitarity limit, 
there are no relevant scale left in the interaction and the dimensionless coefficients bn must be 
universal. This is shown explicitly with a calculation up to the third order in the virial expansion. 
There would be small corrections to these universal results from effective range, higher-partial 



waves, etc., that are neglected in this calculation but are otherwise straightforward to incorporate 
in perturbation. 

The virial coefficient b„ receives contribution from up to and including n-body physics 1121. For 
example, Z?4 would receive contributions from 1-body (non-interacting theory), 2-,3- and 4-body 
physics. b2 is the first coefficient that receives contribution from the interacting theory and it is 
related to the 2-body scattering phase shift 1191 [lOl [Tin : 
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where y= I /a. We use the following notation: 






(3) 



where bn is the n-th virial coefficient from /-body contribution. The thermodynamic pressure Pj^ , 

number density n^ , etc., are also denoted in a similar manner. We also define p(') and n^'^ without 
a subscript to denote the contribution to the pressure and number density from /-body physics to 
all orders in the z expansion. 



Aivjiii j^4.|^ O2 = 1/v^ as a ^ ±00 and it is universal. A plot ofb^ and 
in Fig. [T] near the unitarity limit for lithium ^Li for temperatures around 6^K. 
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FIG. 1: Plot of ^^2 ' the 2-body contribution to the second virial coefficient and its derivative Tdjb^ for 
^Li as a function of the inverse scattering length y= \/a in eV. leVw 5.06 x lO^^A^^ in natural units 
{h= \ = c). Solid, long-dashed, dot-dashed and short-dashed curves correspond to temperatures of 6/7K, 
8.4/iK, 1 1 .76/iK and 16.46/jK respectively (in 40% increments). The second virial coefficient is temperature 
independent and universal as 7= l/a ^ 0*. 



Traditionally the higher virial coefficients are calculated in a cluster expansion lllOll and there is 
no direct simple formula for relating it to many-body interaction. For example, it is not clear how 
the contribution from 3-body bound states associated with Efimov effect 1I12I1 in bosons such as ^He 
near a — > ±0° is taken into account. This is in contrast to an effective field theory calculation where 
the virial coefficients are directly related to many-body scattering and the '^He 3-body bound state 



contribution can be directly taken into account. The general methodology for such an effective 
field theory calculation was developed earlier in Ref. flV\. 

For a dilute system with a clear separation of scales l/\a\ ^ n^'^ <ti ^/MT /{2ii) <ti l/R, ef- 
fective field theory is quite ideal, and provides model independent results in the region of interest. 
A non-relativistic system near the Feshbach resonance in the 5-wave at chemical potential /j (or 
densities n) and temperature T such that y/M/u (or n'/^), ^/MT -^ l/R can be described by the 
Lagrangian density: 
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where the ". . . " represents higher dimensional operators that are suppressed for dilute systems, \\f 
are the spin-^ fermion fields and the o,- matrices act in the spin space. The four-fermion coupling 
g is related to the 2-body scattering length a (see Refs. lll3L ll4ll and references therein): 
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where v is the renormalization scale in dimensional regularization. One can choose to regulate 
using momentum cutoff A and get a ^(A) similar in form to g{\) above [1 J]- The final result is 
independent of the regularization scheme and renormalization scale as it should be. Notice that 
there are no 3 -body 5-wave operator at leading order lll4lll5n . Therefore when we calculate the 
third virial coefficient, it will depend only on the 2-body interaction and will be universal as the 
2-body scattering length a approaches ±0°. 

The organization of the fugacity z expansion in terms of Feynman diagrams was worked out 
earlier in Ref. Illh and applied to bosons. For bosons, there is a 3-body interaction at leading order 
which can lead to a 3-body bound state that dominates the third virial coefficient b^. At the same 
time this prevents b^ from being universal for bosons since the 3-body binding energy is system 
dependent and does not necessarily approach an universal value as a ^ ±00. 

The first two virial coefficients bi and b2 are known. Nevertheless, it is instructive to derive 
them in the effective field theory from Feynman diagrams before calculating the third virial coef- 
ficient bj. Feynman diagrams with a closed particle loop is order z (from the energy integral over 
the Matsubara frequency) and vanish in vacuum as expected. This can be seen by calculating the 
pressure in the free theory which is given by the loop diagram shown in Fig.|2j 
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A closed particle-particle diagram (only possible with interactions) is 0(1) and does not vanish 
in the vacuum. However, a closed loop with "baryon number" equal to 2 i.e. a closed dimer 
propagator is order O(z^). Similarly a closed trimer propagator is O(z^). This is more easily 
demonstrated with the following calculation of the second- virial coefficient. 
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FIG. 2: The pressure P from 1- and 2-body physics. Sohd Unes are fermion fields. Each particle-particle 
loop is order g^/MT ~ a\/MT and summed to all orders to form the dressed dimer propagator, represented 
by the double hnes. 

At large scattering length \a\^yMT /{In) > 1, the dimer propagator is given by an infinite 
geometric sum shown in Fig. El 
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where p = \p\, q = \q\ and pQ is an odd multiple of inT . Some care is needed in evaluating 
^^(POiP) where pQ is to be identified with an "final" loop integral variable in Eqs.|9land[lT] In 
the center of mass frame p = fc — fc = 0, analytically continuing pQ to the center of mass energy 
in vacuum pQ -^ k^ /M — 2/j + zO+, the dimer propagator 'D{k^/M — 2/j, 0) reproduces the correct 
large scattering length 2-body amplitude: 
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The contribution from the dimer propagator to the pressure is 
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The contour integral in T] is in an anti-clockwise sense over all the odd integer multiples of inT . 
Looking at Eqs.|7]and|9l, one can see how the energy integral over T] for a dimer propagator are 
O(z^) because a dimer carries chemical potential 2/j. This is similar to the energy T] integral for a 
closed single particle propagator that carries chemical potential /j and l/[exp(|3ri) + 1] ~ 0(z) for 



T] ~ —/J. P^'^' is O(z^) and an explicit calculation shows: 
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Comparing Eqs.[Tl|2l|6land[T0l we see that the known results up to the second order in the virial 
expansion are exactly reproduced in the effective field theory calculation. 
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FIG. 3: 3-body contribution to the number density n^ 
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The calculation of the third virial coefficient b\ is technically more challenging but it follows 
the same procedure. First, the leading order trimer propagator in the fugacity z expansion is 
calculated, and then the energy integral over imaginary odd multiples of inT is carried out. The 
trimer propagator is 0(1) in the z expansion and carries chemical potential 3/7. Thus an energy 
integration over a closed trimer loop contributes at O(z^). Like the dimer, the trimer propagator 
requires summation of an infinite series of Feynman diagrams for {aly^MT /{In) > 1. However, 
unlike the dimer, the 3-body contributions do not form a geometric series and this makes the 
trimer more challenging. The solution to the trimer problem has been know lll4lll5ll and we will 
not repeat it here. The relevant Feynman diagrams are shown in Fig. |3l It is easier to calculate the 
number density n = d^P instead of the pressure P for the trimer. We find: 
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where the contribution from the trimer propagator, after projecting onto the 5-wave, is defined 
through the integral equations 
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Plot of ^3 and its derivative TdjbT, for ^Li as a function of y = 1/a in eV. The different curves are 
same set of temperatures as in Fig.^and uses the same notations. ^3 shows universal behavior. 



In Fig. 131 the third virial coefficient ^3 = ^3 — Z?, 



^3 + ^3 = ^3 and its derivative Tdrb^ 



are shown for ^Li at temperatures around S/jK. As expected from physical arguments presented 

earlier, the third virial coefficient is universal. This can be made more explicit in Eqs.fTOlandfTTl 

(3) 
for \a\ ^00, and write, for example, the coefficient b^, in terms of dimensionless integrals. We 

find b^ ^ -0.54. 

From this calculation, we see that the virial coefficients up to O(z^) are not un-naturally large 

[numbers of 0(1)] and for z ^ 1, the virial expansion seems consistent. Formally, the expansion 

in z is equivalent to an expansion in diluteness parameter in perturbation ||9lll6|] for nX^ ^ 1: 
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where we used b^ = \. Defining the density n in terms of the Fermi temperature 7p, nh? /2 ~ 
4/{3y/n){Tf/T)^'~ is the small expansion parameter for T ^ Tf. The energy density can be 



obtained from the pressure using standard thermodynamic relations £ 
s = djP is the entropy density. 

As (3 ^ ±00, the energy density can be written as pi 9]: 



-P + p n + T s, where 



3 

2^" 



1 nV n^l^ ( 1 

1 + ..,. ^ -f 



-int 



3 nk^ 

2 2 



25/2 2 

2 + 



M 



4 V8 9V3J 



z.-^) + ^f-v^z.-(^)+4(z.r(2))2 



+ ^int = ^kin+Eint; 
2K 



(14) 



£^j is the contribution from the interacting theory in Eq.|4|in the unitarity limit a -^ ±0°. 

In Fig.|5l we have plotted the contributions from the third virial coefficient to the energy density. 
In the left subplot in Fig.|5J the solid curve shows the contribution up to 0(n?i^) normalized to the 
leading order contribution as a function of Tp/T . The dashed curve is the O(n^X^) contribution 
which is a smaller correction for Tp/T 5 0.49. This is shown more clearly in the right subplot in 
Fig. |5J the convergence (error) plot HlTll . It is deduced that the range of convergence is Tf/T < 
0.49. 

The ratio of the interaction energy 8int to the contribution from the free Lagrangian density 
Ekin is shown in Fig. |6lfor a range of Tf/T where the virial expansion converges at this order of 
the calculation. For the experiments on ^Li, at temperatures around T ~ 3Tf we expect the virial 
expansion to be valid. At this temperatures, the effect of the third virial coefficient should be 
detectable in experimental data for dilute systems satisfying the constraint n^'^ <^ -^MT j{2'K) <^ 
\/R. 
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FIG. 5: Left subplot: the second (solid curve) and third order (dashed curve) energy density £°° in the n'k} /2 
expansion in Eq.Q4]as a fraction of the leading order result 'iTn/2 vs. Tp/T . Right subplot: Convergence 
(eiTor) plots ilTT . Absolute value of the ratio of only the Qinkr') to the leading order contribution (solid 
curve) and only the 0{n^X^) to the leading order contiibution (dashed curve) to the energy density as a 
function of Tf/T . For small values of Tf/T , the 0{n^X^) effects are smaller than 0{n}?). These effects 
are equal at around T ~ 2Tf and the perturbation breaks down. Note the 0(«^A,^) effects equal the leading 
order result (ratio equals 1) at a smaller temperature T ~ \Tf, at fixed Tf. The O(nX^) results are smaller 
than the leading order results to even smaller temperatures (not shown). 

In conclusion, we considered a spin-^ dilute fermionic system near the unitarity limit \a\ —>■ 0°. 
At finite temperature the thermodynamic pressure was calculated in a model independent way 
using effective field theory up to the third order in the virial expansion. At this order, universality 
was demonstrated and a range of temperature T >3Tf for a fixed density (conversely a range of 
densities for a fixed temperature) was identified for a self-consistent application of the results. 
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FIG. 6: Ratio sj^j/Sjt,,, of interaction to kinetic energy as a function of Tp/T in the range of convergence 
Tf/T < 0.49 deduced in Fig. |5] Solid curve: leading order [0(?iX^)] result, dashed curve: next-to-leading 
order [0{n^X^)] result. 

The virial expansion is an useful tool for dilute systems where the microscopic physics is poorly 
known or difficult to calculate. The universal forms for the equation of state derived here are 
applicable to a wide class of physical problems in nuclear and atomic physics. These results 
are relevant for understanding neutrino physics in hot dilute nuclear matter n ~ 10^^ fm-' 1^, 
and the same equations for energy density, etc., are testable in atomic physics experiments. At 



lower temperatures, numerical lattice calculations are more appropriate II18LI19LI20I1 . However, at 
intermediate temperatures, our results could be used as checks for the lattice calculations. 
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